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Following a previous detailed study of unpolarized generalized parton distri- 
bution functions in the non-singlet sector, and at zero values of the skewness 
variable, (, we propose a physically motivated paramctrization that is valid 
at £ 7^ 0. Our method makes use of information from the nucleon form factor 
data, from deep inelastuc scattering parton distribution functions, and from 
lattice results on the Mellin moments of generalized parton distributions. It 
provides, therefore, a step towards a model independent extraction of gener- 
alized distributions from the data, alternative to the mathematical ansatz of 
double distributions. Comparisons with recent experimental data on the proton 
are shown. 



1. Introduction 

Generalized Parton Distributions (GPDs) parametrize the soft contribu- 
tions in a variety of hard exclusive processes, from Deeply Virtual Compton 
Scattering (DVCS) to hard exclusive meson production (see 1,2 for reviews). 
The feasibility of DVCS-type experiments using a deep inelastic probe with 
virtuality, Q 2 , while detecting a momentum transfer, t, between the ini- 
tial and final proton, allows one to address a vast, previously inaccessible 
phenomenology. In particular, one can access parton densities in impact 
parameter space, 3 and envisage extracting the orbital angular momentum 
of partons in both nucleons and nuclei. 4 

At present, a central issue is the definition of a quantitative, reliable 
approach beyond the construction of GPDs from specific models and/or 
particular limiting cases, that can incorporate new incoming experimen- 
tal data in a variety of ranges of Q 2 , and of the longitudinal (along the 
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lightcone) and transverse components of the four-momentum transfer be- 
tween the incoming and outgoing protons, ( and Aj_, respectively (see also 
discussion in 5 ). The matching between measured quantities and Perturba- 
tive QCD (PQCD) based predictions for DVCS should proceed, owing to s 
pecific factorization theorems, similarly to the extraction of Parton Distri- 
bution Functions (PDFs) from deep inelastic scattering. A few important 
caveats are however present since GPDs describe amplitudes and are there- 
fore more elusive observables. The comparison between experiment and 
the formulation of parametrizations necessarily encompasses other strate- 
gies using additional constraints, since a direct comparison with the data 
cannot be performed. Experiments providing sufficiently accurate data to 
constrain the shape of GPDs have just begun. 6,7 

In Refs. 8,9 we proposed a strategy using a combination of experimental 
data on nucleon form factors, PDFs, and lattice calculations of Mellin mo- 
ments with n > 1. The latter, parametrized in terms of Generalized Form 
Factors (GFFs), were calculated by both the QCDSF 10 and LHPC 11 col- 
laborations for both the unpolarized and polarized cases up to n — 3. At 
n > 2, due to the polynomiality property, 1 ' 2 the Mellin moments become 
dependent on the skewness, £. 

In this contribution we report on the approach used in; 9 we show our 
results for both proton and neutron DVCS, and we finally present a pre- 
liminay result on the angular momentum sum rule in the deuteron. 

2. Generalized Parton Distributions from Lattice Moments 

In this contribution we concentrate on the unpolarized scattering GPDs, 
H, and E, from the vector (7^) and tensor (a^ u ) interactions, respec- 
tively. We adopt the following set of kinematical variables: ((, X, t), where 
( = Q 2 /2(Pq) is the longitudinal momentum transfer between the ini- 
tial and final protons (( « Xsj in the asymptotic limit, with Bjorken 
%Bj = Q 2 /2Mv), X = (kq)/(Pq) is the momentum fraction relative to the 
initial proton carried by the struck parton, t = — A 2 , is the four-momentum 
transfer squared. X is not directly observable, it appears in the amplitude 
as an integration variable. 1,2 The need to deal with a more complicated 
phase space, in addition to the fact that DVCS interferes coherently with 
the Bethe-Heitler (BH) process, are in essence the reasons why it is more 
challenging to extract GPDs from experiment, wherefore guidance from 
phenomenologically motivated parametrizations becomes important. 
We first present a paramctrization of H and E in the flavor Non Singlet 
(NS) sector, valid in the X > £ region, obtained by extending our previous 
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zero skewness form, 8 through proper kinematical shifts: 



H(X,C,t) = G x Mx (X,C,t)R(X,C,t) 



(1) 



(a similar form is obtained for E(X, £, t)) 1 where R(X, £, t) is a Regge mo- 
tivated term describing the low X and t behaviors, while G\ Ix (X, £, t), was 
obtained within a spectator model. 

In order to model the X < ( region, we observe that the higher moments 
of GPDs give ^-dependent constraints, in addition to the ones from the 
nucleon form factors. The n — 1,2,3 moments of the NS combinations: 
H u-d = H u_ H d^ and E u-d = E u_ E d are avai i a ble from lattice QCD. 10 ' 11 

They can be written in terms of the isovector components as: 



where the l.h.s. quantities are obtained from the lattice moments calcula- 
tions, whereas (H^j^jn and (E^^^n are amenable to chiral extrapola- 
tions. We used lattice calculations for the unpolarized GFFs obtained by 
the QCDSF collaboration using two flavors of 0(a)-improved dynamical 
fcrmions for several values of t in the interval < t < 5 GcV 2 , and cov- 
ering a range of pion mass values, m w > 500 MeV 2 . Similarly to previous 
evaluations 11 the GFFs for both H and E, display a dipole type behavior 
for all three n values, the value of the dipole mass increasing with n. We 
performed an extrapolation by extending to the n = 2, 3 moments a simple 
ansatz proposed in 12 for the nucleon form factors that: i) uses the con- 
nection between the dipole mass and the nucleons radius; ii) introduces a 
modification of the non analytic terms in the standard chiral extrapolation 
that suppresses the contribution of chiral loops at large m T . Despite its 
simplicity, the ansatz seems to reproduce both the lattice results trend at 
large m T while satisfying the main physical criteria i) and ii). Our results 
for the dipole mass at n = 2 are shown in Fig.l. 

3. Reconstruction from Bernstein Polynomials 

Similarly to the PDFs case, 13 with a finite number of moments in hand, one 
can use reconstruction methods attaining weighted averages of the GPDs, 
around average ranges of X. The weights are provided by the complete set 
of Bernstein polynomials. The Bernstein polynomials are ideal for repro- 
ducing the deep inelastic structure functions in that they are zero at the 
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Fig. 1. The dipole masses squared for n = 2, for the isovector magnetic (lower panel) 
and electric (upper panel) contributions obtained by performing fits to the lattice results 
of. 10 The value at the physical pion mass obtained from our fit is also shown (star), 
(adapted from 9 ) 



endpoints, they are normalized to one, and they are peaked in different 
regions within the interval xbj £ [0, f]. Because of the latter property the 
Bernstein polynomials allow one to emphasize the behavior of the structure 
function at given specific regions of XBj, while suppressing the others. It 
was found that n > 8 moments were necessary to give a fully quantitative 
description of the behavior of F2(xbj,Q 2 )- At present, only n — 3 Mellin 
moments are available from lattice QCD calculations, therefore one cannot 
reach a similar level of accuracy as in the inclusive case. 

In Fig. 2 we show H u ~ d reconstructed using the available lattice mo- 
ments. We performed the procedure in the X < ( region only using: 

C 

H k ,n(C,t) = J H(X,C,t)b kin (X,QdX k = 0,...n, (4) 
o 

C 

where: b k ,n{ X '0 = X k {( - X) n - k / J X k (Q - X) n - k dX, and we used 

o 

subtracted moments, defined as: 

l 

( H n)x<c= H n- [ H\X,C,t)X n dX, (5) 

c 

where H n are the Mellin moments, and H 1 (X,(,t) was obtained from 
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Eq.(l). For n 



Hi2((Xi2) 



H22{(,X22) 

where X i — 0.25, X 02 = 0.5, X 3 = 0.75, and A w , A 2 o, A 30 , A 32 are the 
GFFs from Ref. 10 





Fig. 2. Comparison of H u ~ d for different values of ( = 0.18, 0.25, 0.36, 0.53, and -t = 
tmin = 0.035, 0.073, 0.18, 0.53 GeV 2 calculated using the procedure described in the text 
(adapted from 9 ). 



= 2, k = 0, 1, 2, the reconstruction procedure yields: 
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(6) 



4. Comparison with Experiment 

The unpolarized Compton form factors were recently extracted from DVCS 
on proton and deuteron targets at Q 2 « 2 GeV 2 , C = x Bj = 0.36, for several 
values of t in the range 0.15 < t < 0.35 GeV 2 . 6,7 In Fig. 3 and we show the 
results of our parametrization at C — X for the imaginary part of the 
Bethe-Heitler BH-DVCS interference term at leading order, for proton and 
neutron, respectively. The proton results are compared with recent data 
from Jefferson Lab, 6 at xbj = C = 0.36, and Q 2 w 2 GeV 2 , while the 
neutron ones are presented at the kinematics of the forthcoming analysis 
from Ref. 7 . 
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Fig. 3. The imaginary part of the Bethe-Heitler BH-DVCS interference term at leading 
order, C I (F), (see 1 ), (upper panel) Proton target. Experimental data from Ref. 6 The 
full curve shows our prediction, including the theoretical error, evolved to Q 2 = 2 GeV 2 . 
The dot-dashed curve, obtained at Q 2 = 4 GeV 2 , shows that the effect of evolution is 
relatively small. The other curves represent variations of the parametrization obtained 
respectively by disregarding the contribution from E, by disregarding the kincmatical 
shift in the Regge term described in the text, and by using E from Set II of Ref. 8 
(lower panel) Neutron target. The two curves show our results using variations of the 
parametrization for E using Set I and Set II of Ref. 8 The curves are in agreement with 
the experimental data from Ref. 7 (adapted from Ref. 9 ) 



5. Angular Momentum Sum Rule 

We present results derived for the angular momentum for a spin 1 system 
in terms of second moments of it's generalized parton distributions (a full 
derivation will be shown in 14 ). 

The second moment of the parton distributions is relevant to the spin 
structure of the nucleon. To see this, we first write down the angular mo- 
mentum operator in QCD as the sum of quark and gluon contributions, 15 



where 



Jqcd — J q + J g 



d 3 x x x T„ 



(7) 



d 3 a 



2-j — * 

ip^—ip + x (-iD)'ip 



Jo 



d 3 x xx(E x B) 



(8) 



Here T q and E x B are the quark and gluon momentum densities, respec- 
tively. S is the Dirac spin-matrix and D = d + igA is the covariant deriva- 
tive. By an analogy with the magnetic moment, one can get the separate 
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quark and gluon contributions to the nucleon spin if the form factors of the 
momentum density, or equivalently the energy-momentum tensor of QCD, 
are known at zero momentum transfer. Using Lorentz covariance and other 
symmetry principles, one can write down six form- factors separately for 
quark and gluon parts of the energy-momentum tensor, 

= -\ \P"P V ] (e'*e)G li2 (i) - \ [P»P V ] {tP) ^* P) G 2a (t) 

- I [A"A"- 5 '«'A 2 ] ( c '* c )G 3 ,2(t) 

_ 1 [A-A'-^A 2 ] {tP) ^ P) G^(t) 

+ i [(e'*"(eP) + e"(e'*P)) P v + /i <-> i/] G 5 , 2 (i) 
+ [(e'* M (eP) - e^(e'*P)) A^ + /i «-> 1/ 

+ 2g^(eP)(e'*P) - (e"^e' y + e'^e") A 2 ] G 6;2 (*) (9) 

whereP M = p^ + p^' , A p = — and e is the deuteron polarization. 
Taking /i = 0, in the Brcit frame, the matrix element of J q , gives 

(p'\ J d 3 x(x x f% g ) z \p) = G 5 , 2 (0) | d 3 x p° (10) 

G 5;2 (0) = 2J« (11) 

Therefore, the generalized form factor, G5 >2 (0), of the energy momentum 
tensor of spin 1 system provides the fraction of the spin is carried by quarks 
and gluons. 

6. Conclusions 

In conclusion, we provided a fully quantitative parametrization of the GPDs 
in the non-singlet sector, valid in the region of Jefferson Lab experiments 
using both proton and deuteron targets. 6 ' 7 We also presented a new result 
for the angular momentum sum rule in the deuteron. Differently from model 
calculations, and for the first time to our knowledge, our parametrization 
makes use of experimental data in combination with lattice results. Given 
the paucity of direct experimental measurements of GPDs, and the uncer- 
tainties related to the partonic interpretation of GPDs in the ERBL region, 
our goal is to provide more stringent, model independent predictions that 
will be useful both for model builders, in order to understand the dynamics 
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of GPDs, and for the planning of future hard exclusive scattering experi- 
ments. 
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